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RESUMEN 

El resumen sera traducido al espanol por los editores. Using numerical simulations and the Sheth- 
Tormen approximation we study the mass function of dark matter halos in voids. We find that the void mass 
function is significantly lower and its shape is different than that of the field halos. We predict that in the 
standard LCDM model a void with radius 10ft.~ 1 Mpc should have 50 halos with circular velocity v c > 50 km/s 
and 600 halos with v c > 20 km/s. 

ABSTRACT 

Using numerical simulations and the Sheth-Tormen approximation we study the mass function of dark matter 
halos in voids. We find that the void mass function is significantly lower and its shape is different than that 
of the field halos. We predict that in the standard LCDM model a void with radius 10/i _1 Mpc should have 50 
halos with circular velocity v c > 50 km/s and 600 halos with v c > 20 km/s. 
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According to the hierarchical structure formation 
paradigm halos form from small initial density fluc- 
tuations by accretion and merging. It is expected 
that inside the dark matter (DM) halos gas cools and 
galaxies are formed. Thus, we expect that the dis- 
tribution of DM halos and their structure are closely 
related to the observed galaxy distribution. It is also 
known (e.g., Gottlober et al. 2001) that the merging 
history of the DM halos depends on environment. 
Here we study DM halos in the low density environ- 
ment of voids. 

Cosmological simulations predict many more 
small DM halos than the observed number of satel- 
lites around the Milky Way and Andromeda galax- 
ies. Do we have the same problem for dwarf galaxies 
in voids? One naively expects a large number be- 
cause the Press-Schechter mass function steeply rises 
with declining mass. In contrast, it seems that ob- 
servations are failing to find a substantial number 
of dwarf galaxies inside voids (e.g. Popescu et al. 
1997). However, the situation is complicated because 
it is very difficult to detect dwarf galaxies, many of 
which are expected to have low surface brightness. 

We use N-body simulations to make accurate pre- 
dictions for the expected number of dwarf halos in 
voids. We have performed a series of high resolu- 
tion simulations using the ART code (Kravtsov et 
al. 1997) and the friends-of-friends (FOF) algorithm 
to identify halos. We investigate a spatially flat 



cold dark matter model with a cosmological constant 
(ACDM with Q M = 0.3, Q A = 0.7, ct 8 = 0.9, and h 
0.7 We use a cube of 80/i _1 Mpc side length, which 
is sufficiently large to study the formation of large 
voids. Because we are interested in the formation 
of small structure elements inside the voids, we need 
high mass resolution. Therefore, we perform the sim- 
ulation in two steps. First, we run a low mass reso- 
lution simulation (128 3 particles with mass m pa rt = 
2.0 x I0 w h- 1 M Q ) and identify 8387 galactic-size ha- 
los with masses > 2.0 x 10 1 1 /i~ 1 M Q . This corre- 
sponds to a mean distance of about 4/i _1 Mpc be- 
tween halos. We identify voids in the distribu- 
tion of the halos. The largest void has a diame- 
ter of 24/i" 1 Mpc. Second, we re-simulate the voids 
with a formal mass resolution of 1024 3 particles, i.e. 
T^part = 4.0 x 10 7 /i^ 1 Af Q . Thus, we resolve in the 
voids objects with masses larger than 10 9 /i _1 Mq. 
By construction, the voids do not have halos with 
masses larger than M h = 2.0 x 1O 11 /i~ 1 M . 

We resimulate five voids with radii i? V oid = 
11.6, 10.8, 9.4, 9.1, 9.1 ft, _1 Mpc. The underdensity at 
the centers of these voids is 0.1 - 0.2 of the mean 
matter density CIm- I n Fig. 1 we compare the mass 
function of DM halos in voids with the mass function 
in the whole box. The thin solid line is the analyti- 
cal prediction based on the Sheth and Tormen (1999) 
correction to the Press-Schechter mass function 
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Fig. 1. Mass function of DM halos. The thick full 
line shows the cumulative mass function N(> M) in 
the whole 80h~ Mpc simulation. The thin full curve is 
the Sheth-Tormen approximation. The mass functions 
in voids are shown by thick dashed curves. The thin 
dashed curve is our Sheth-Tormen-type prediction of the 
void mass function. The mass functions of the voids are 
different because of the difference in the void underden- 
sity. Note that the normalization and the shape of the 
void functions are different from that of the field popu- 
lation. 

where A = 0.322, p = 0.3, a = 0.707, g h is the back- 
ground density, a is the rms density fluctuation on 
scale R corresponding to mass M (M — Ang^R 3 /3) 
and S c = 1.676 for our model. As expected, within 
the whole mass range covered by the numerical mass 
function it agrees well with the analytical prediction. 
Small deviations at M > 10 14 /i _1 M© are due to low 
number of massive objects. 

The simulated voids have somewhat different 
mean (under)dcnsities. This results in different mass 
functions. The number-density of halos in voids is 
about an order of magnitude smaller than in the 
whole box as expected due to lower mean density 
in voids. However, the shape of the mass function 
is also different. An analytic approximation to mass 
function in voids, shown as a thin dashed line in Fig. 
1, was obtained by applying the formalism of Shcth 
and Tormen (2002) for the constrained mass func- 
tions 
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The parameter So is the linear underdensity of the 
void corresponding to the actual nonlinear under- 
density 5 and is calculated from the spherical top- 
hat model, while <jq is the rms fluctuation at scale 
Ro = (1 + Sy^Rvoid- The parameters a = 0.615, 
(3 = 0.485 come from the ellipsoidal dynamics and 
we have chosen a = 0.5 which fits the simulated 
mass functions better than the value of a = 0.707 
advertised by Sheth and Tormen. The result shown 
in Fig. 1 was obtained for f2 vo id — 0.04 and i? vo id — 
10/i _1 Mpc. In an upcoming paper (Gottlober et al. 
2002, in preparation) we discuss the spatial distribu- 
tion of halos found in numerical simulations of voids 
and their mass function, and we present a straight- 
forward analytical method to predict the mass func- 
tion. 

In a typical void of radius 10/i _1 Mpc we found 
about 600 halos with circular velocities larger than 
20 km/s and about 50 halos with circular velocities 
larger than 50 km/s. Testing these predictions could 
be one of the interesting observational tasks of the 
ongoing SLOAN Digital Sky Survey. 
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